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Abstract. We briefly overview our newly developed radiation transport module for MHD sim- 
ulations and two actual applications. The method combines the advantage of the speed of the 
Flux-Limited Diffusion approximation and the high accuracy obtained in ray-tracing methods. 
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1. Radiative MHD 

Aim of the development of the radiation transport module described here is to achieve 
a fast method for approximative frequency-dependent radiation transport for (magneto-) 
hydrodynamical simulations in studies of the environment around a single object source. 

The radiation transport module provides a Flux-Limited Diffusion (hereafter FLD) 
solver in cartesian, cylindrical and spherical coordinates with a potentially non-equidistant 
grid spacing in the first dimension, e.g. logarithmic in the radial coordinate. Addition- 
ally to the FLD approximation the module consists of a first order ray-tracing technique 
along the first dimension to account for single source irradiation (e.g. stellar heating of 
a proto-planetary disk or radiative feedback in massive star formation). 

The local radiation field is split into two parts, originated from an externally irradiated 
(frequency-dependent) flux Fiy) and a diffuse radiation energy density Ep, which are in 
equlibrium with the radiation from the dust grains: 

aT'^ = Er + 3—- / k{v) F{v) dv, (1.1) 
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where a is the radiation constant, c is the speed of light, k(i') and Kp(T) represent the 
frequency-dependent and the Planck mean opacity for a given temperature T respectively. 
Gas and dust temperatures are assumed to be the same. The irradiation is treated as 
an instantaneous source of additional energy, calculated via a first order ray-tracing as a 
function of the optical depth T^(r) and distance r from the central star: 

F{iy,r) = Kii^) (^) e-^'-W (1.2) 



The evolution of the radiation energy is described by a Flux-Limited Diffusion equation 
dtEp^ = -fc (^V • (^D^Er + J F{v, r)dv^ - (1.3) 

with fc = {cyp/iaT^ + 1)^^, D = Xc/npp and = — pV • u-|-additional source terms 
from hydrodynamics. The flux-limiter A is chosen according to |Levermore fc Pomraning| 
[]( |1981[ ), Cy is the specie heat capacity, p, p and v the gas density, thermal pressure and 



dynamical velocity and kr specifies the Rosseland mean opacity. Scattering is neglected. 
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2. Application I: Massive Star Formation 

Due to the short Kelvin-Helmholtz contraction timescale of a massive star ( Shu et al 
[] (1987)) the accretion process onto such a star is described by the interaction of the 



gravitationally forced inflow of matter with the radiative force escaping from the newly 
born star. The conservation of angular momentum leads to the formation of an accretion 
disk as well as polar cavities. At present we're studying this collapse scenario with respect 
to the effect of dimensionality (from 1 to 3D) and different applied physics (isothermal 
and adiabatic test runs, realistic cooling and frequency-averaged as well as frequency- 
dependent radiative feedback). 



3. Application II: Accretion Disks 

We perform accretion disk simulations with a hydrodynamical stable stratified disk 
and an initial toroidal magnetic field. The initial plasma beta (3 is 



model (f = 0.05 
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constant everywhere in the disk at a value of 25, the runs by Fromang & Nelson 
(hereafter FN). The convergence tests have a radial extension from 4 to 6 AU, 
scale height as vertical and 7r/3 as azimuthal extension. We choose perodic boundaries 
in vertical and azimuthal direction. Small velocity perturbations in radial and vertical 
direction about 10~'*co initiate the non-linear MRI evolution. With the used high-order 



MHD Rieniann solver (Mignone, et al. (20071) we get a converged alpha value of about 



10^^ presented in figure [ij In opposite, FN got an increasing alpha value from low 
a = 10"*^ to high a = 5 * 10~^ (yet still not converged) with increasing resolution for a 
non-Riemann solver. In the future we will use the radiative module to calculate proper 
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Figure 1. Left: Azimuthally averaged logarithmic plasma beta in the turbulent phase. 
Right: Evolution of Maxwell alpha at different resolution, converging against each other. 

resistive terms rj{T, p) and to handle correctly the resulting heating and cooling processes 
of MHD turbulence and the behaviour for the alpha stresses. 
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